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(57) Abstract: Load-Flow computations are performed in 
real-time operation/control and in on-line/off-line studies 
of electrical power systems. Three Load-Flow computation 
methods of the present invention are the best versions of many 
simple variants with almost similar performance. The Super 
Super Decoupled Loadflow (SSDL-YY) method and its many 
variants are characterized in limiting rotation angle applied to 
nodal real and reactive power mismatches to the maximum of 
-48 degrees instead of -36 degrees, replacing network shunt 
element -2bpCos* p by {2(QSH p Cos* p - PSHpSin^pVV, 1 ) or 
by [-bpCosfrp + (QSH p Cos* p - PSHpSin* p y\\ 2 ) and using the 
dividing term V 2 instead of V in the modified nodal real power 
residues [RP] in the prior art Fast Super Decoupled Loadflow 
(FSDL) method. The other two Super Super Decoupled 
Loadflow: BGX' version (SSDl^BGX^and X'G^X* version 
(SSDL-X'GprX') are characterized in the use of simultaneous 
(IV, 18) iteration scheme thereby reducing the mismatch 
computation once compared to two mismatch computations 
in the prior art method employing successive (18,1V) iteration 
scheme. The invented methods are also characterized in the 
different definition of gain matrices as detailed in the steps of 
algorithm -2, algorithm-3 and algorithm-4 of carrying out of 
the inventions, (steps-cc, -dd, and ff in Fig.2; steps-ccc, -ddd, 
-ffi". -888 am* -hhh m R 8-3; steps-dddd and -ABET in Fig. 4) 
leading to some speed-up of the invented methods. 
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SYSTEM OF SUPER SUPER DECOUPLED LOADFLOW COMPUTATION FOR 

ELECTRICAL POWER SYSTEM 

TECHNICAL FIELD 

The present invention relates to methods of loadflow computation in power flow control 
and voltage control for an electrical power system. 

» 

BACKGROXJND ART AND MOTIVATION 

Utility/industrial power networks are composed of many power plants/generators 
interconnected through transmission/distribution lines to other loads and motors. Each of 
these components of the power netwoik is protected against unhealthy (faulty, over/under 
voltage, over loaded) potentially damaging operating conditions. Such protection is 
automatic and operates without the consent of power network operator, and takes an 
unhealthy component out of service disconnecting it from the network. The time domain of 
operation of the protection is of the order of milliseconds. 

The purpose of a utility/industrial power netwoik is to meet the electricity demands of its 
various consumers 24-hours a day, 7-days a week while maintaining the quality of 
electricity supply. The quality of electricity supply means the consumer demands be met at 
specified (say + or - 5% tolerance) voltage and frequency levels without over loaded, 
under/over voltage operation of any of the power netwoik components. The operation of a 
power netwoik is different at different times due to changing consumer demands and/or 
development of any faulty/contingency situation. In other words healthy operating power 
network is constantly subjected to small or large disturbances. These disturbances could be 
operator initiated, or initiated by security control functions and various optimization 
functions such as economic operation, minimization of losses etc., or caused by a 
fault/contingency incident 



SUBSTITUTE SHEET (RULE 26) 



WO 2004/023622 PCT/CA2003/001312 

For example, a power network is operating healthy and meeting quality electricity needs of 
its consumers. A fault occurs on a line or a transformer or a generator which faulty 
component gets isolated from the rest of the healthy network by virtue of the automatic 
operation of its protection. Such a disturbance would cause a change in the patten of 
power flows in the network, which can cause over loading of one or more of the other 
components and/or over/under voltage at one or more nodes in the rest of the network. This 
in turn can isolate one or more other components out of service by virtue of the operation 
of associated protection, which disturbance can trigger chain reaction disintegrating the 
power network. 

Therefore, the most basic and integral part of all other functions (e.g. optimizations) in 
power network operation and control is security control. Security control means controlling 
power flows so that no component of the netwoik is over loaded and controlling voltages 
such that there is no over voltage or under voltage at any of the nodes in the network 
following a disturbance small or large. Security control functions (overload alleviation and 
over/under voltage alleviation) can be realized through one or combination of more 
controls in the network These involve control of power flow over tie line connecting other 
utility netwoik, turbine steam/water input control to control real power generated by each 
generator, load shedding function curtails load demands of consumers, excitation controls 
reactive power generated by individual generator which essentially controls generator 
terminal voltage, transformer taps control connected node voltage, switching in/out in 
capacitor/reactor banks controls reactive power at the connected node. Such overload and 
under/over voltage alleviation functions produce control amount changes in step-60 of 
Fig. 5. These control amount changes could be even optimized in case of simulation of all 
possible imaginable disturbances (outage of a line, loss of generation etc.) for corrective 
action stored and made readily available for acting upon in case the simulated disturbance 
actually occurs in the power network. In fact simulation of all possible imaginable 
disturbances is the modern practice because corrective actions need be taken before the 
operation of individual protection of unhealthy component 

2 
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Control of an electrical power system (Power-flow control, voltage control etc.) is 
performed according to the process flow diagram of fig. 5. The various numbered steps in 
fig. 5 are explained in the following. 

Step- 10: On-line readings of various real-time power flows, voltages, circuit breaker status 

(open/close) etc. are obtained 
Step-20: A control amount (Le. change in power injections, voltages etc.) is initially 

established and proposed 
Step-30: Various power flows, voltage magnitudes and angles, reactive power generations 

by generators and turns ratios of transformers in the power system are determined 

by performing loadflow computation, which incorporates established/proposed/set 

control adjustments 

Step-40: The results of Loadflow computation of step 30 are evaluated for any over loaded 

transmission lines and over/under voltages at different nodes in the power system 
Step-50: If the system state is good (no over loaded lines and no over/under voltages), the 

process branches to step 70, otherwise to 60 
Step-60: Changes the control amount initially set in step-20 or later set in the previous 

process cycle step-60 and returns to step-30 
Step-70: Actually implements the control amount correction to obtain 

secure/optimum/correct/acceptable operation of power system 

It is obvious that Loadflow computation is performed many times in real-time operation 
and control environment and, therefore, high-speed (efficient) Loadflow computation is 
necessary to provide corrective control in the changing power system conditions including 
an outage or failure. Moreover, the loadflow computations must be highly reliable to 
yield converged solution under wide range of system operating conditions and 
network parameters. Failure to yield converged loaflow solution creates blind spot as to 
what exactly could be happening in the network leading to potentially damaging 
operational and control decisions/actions in capital-intensive power utilities. 
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The embodiment of the present invention, the most efficient and reliable loadflow 
computations, as described in the above steps and in Fig. 5 is very general and elaborate. 
The control of voltage magnitude within reactive power generation capabilities of electrical 
machines (generators, synchronous motors, capacitor/inductor banks) and within operating 
ranges of transformer taps is normally integral part of Loadflow computation as described 
in "LTC Transformers and MVAR violations in the Fast Decoupled Loadflow, IEEE PAS- 
101, No.9, PP. 3328-3332." If under/over voltage still exists in the results of Loadflow 
computations, other control actions are taken in step-60 in the above and in Fig. 5. For 
example, under voltage can be alleviated by shedding some of the load connected. 

However, the simplest embodiment of the efficient and reliable system and method of 
loadflow computations is where only voltage magnitudes are controlled by controlling 
reactive power generation of generators and motors, switching in/out in capacitor/inductor 
banks and transformer taps. Of course, such control is possible only within reactive power 
capabilities of machines and capacitor/reactor banks, and within operating ranges of 
transformer taps. This is the case of a network in which the. real power assignments have 
already been fixed and in which steps-50 and -60 in the above and in Fig. 5 are absent Once 
loadflow computations are finished, the Loadflow solution includes indications of reactive 
power generation at generator nodes and at the nodes of the capacitor/inductor banks, and 
indications of transformer tap settings. Based on the known reactive power capability 
characteristics of the individual machines, the determined reactive power values are used to 
adjust the excitation current to each machine, or at least each machine presently under 
reactive power, or VAR, control, to establish the desired reactive power set points. The 
transformer taps are set in accordance with the tap setting indications produced by the 
Loadflow computation system. 

This procedure can be employed either on-line or off-line. In off-line operation, the user 
can simulate and experiment with various sets of operating conditions and determine 
reactive power generation and transformer tap settings requirements. A general-purpose 
computer can implement the entire system For on-line operation, the loadflow computation 
system is provided with data identifying the current real and reactive power assignments 
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and transformer transformation ratios, the present status of all switches and circuit breakers 
in the network and machine characteristic curves in steps- 10 and -20 in the above and in 
Fig. 5, and blocks 10, 12, 14, 20, 30, 40, 42, 44, 50, 52 in Fig 6. Based on this information, 
a model of the system based on gain matrices of any of the invented or prior art Loadflow 
computation methods provide the values for the corresponding node voltages, reactive 
power set points for each machine and the transformation ratio and tap changer position for 
each transformer. 

■ 

The present invention relates to control of utility/industrial power networks of the types 
including plurality of power plants/generators and one or more motors/loads, and 
connected to other external utility. In the utility/industrial systems of this type it is the usual 
practice to adjust the real and reactive power produced by each generator and each of the 
other sources (synchronous condensers, capacitor/inductor banks) in order to optimize the 
real and reactive power generation assignments of the system. Healthy (secure) operation 

■ 

of the network can be shifted to optimized operation through corrective control 
(disturbance) produced by optimization functions without violation of security constraints. 
This is referred to as security constrained optimization of operation- Such an optimization 
is described in the United States Patent Number. 5,081,591 dated Jan. 13, 1992 
(Optimizing Reactive Power Distribution in an Industrial Power Network) where the 
present invention can be embodied by replacing the block nos. 56 and 66 by a block of 
constant matrices [Y9] and [YVJ, and replacing blocks of ''Exercise Newton-Raphson 
Algorithm" by blocks of "Exercise Fast Super Decoupled Algorithm" or "Execise Super 
Super Decoupled Algorithm" in place of blocks 58 and 68. 

DISCLOSURE OF THE INVENTION 

This invention relates to steady-state power network computation referred to as Loadflow 
or Power-Flow. Loadflow computations are performed as a step in real-time 
operation/control and in on-line/off-line studies of Electrical Power Systems. The present 
invention involves three-methods. These 3-methods are the best versions of many simple 
variants with almost similar performance. Simple variants include any possible hybrid 
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combination of these 3 -methods and u asymmetrical definitions of [Y0] in SSDL-method. 



Among these 3-methods, their variants and all other known methods, Super Super 
Decoupled Loadflow (SSDL-YY) is the simplest, easiest to implement and overall best 
in performance (reliability of convergence and efficiency of computations). 

BRIEF DESCRIPTION OF DRAWINGS 

i 

Fig. 1 is a flow-chart of the prior art Fast Super Decoupled Loadflow computation method 
Fig. 2 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow 

computation method of version SSDL-YY 
Fig. 3 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow: 

SSDL-BGX\ -BGY, and -BGX versions 
Fig. 4 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow: 

SSDL-X'GpvX\ SSDL-YGpvY, and SSDL- XGpyX versions 
Fig. 5 is a flow-chart of the overall controlling method for an electrical power system 

involving Loadflow computation as a step which can be executed using one of the 

Loadflow computation methods of Figs. 1, 2, 3, 4, other variations described, their 

* * 

hybrid combination and/or their simple variants 
Fig. 6 is a flow-chart of the simple special case of voltage control in overall controlling 
method of Fig. 5 for an electrical power system 



Symbols 



The prior art and inventions will now be described using the following symbols: 



Y pq~Gpq + jB pq 
V p = ep+jf p =VpZ0p 
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P P + jQp 

APp+jAQp 

RPp+jRQp 
PSHp+jQSHp 



k 



m 



net nodal injected power, calculated 

nodal power residue (mismatch) 

modified nodal power residue 

net nodal injected power, scheduled 

rotation angle 

number of PQ-nodes 

number of PV-nodes 



n=m+k+l 
q>p 



[ 1 

LRA 



total number of nodes 

q is the node adjacent to node-p excluding the case of q=p 
indicates enclosed variable symbol to be a vector or a matrix 
Limiting Rotation Angle 



PQ-node 
PV-node 



: load-node (Real-Power-P and Reactive-Power-Q are specified) 



: generator-node (Real-Power-P and Voltage-Magnitude- V are specified) 



Decoupled Loadflow 

A class of decoupled Loadflow methods involves a system of equations for the separate 
calculation of voltage angle and voltage magnitude corrections. Each decoupled method 
comprises a system of equations (1) and (2) differing in the definition of elements of [RP], 
[RQ], and [Y9J and [TV]. 



In this scheme (1) and (2) are solved alternately with intermediate updating. Each iteration 



[RP] = [Y8] [A9] 



(1) 



IRQ] = IYV] [AV] 



(2) 



Successive (18, IV) Iteration Scheme 
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[A9] = [Y9] - 1 [RP] 



(3) 



[e] = [e] + m 



(4) 



[AV] = [YV] 1 [RQ] 



(5) 



[V] = [V] + [AV] 



(6) 



The scheme involves solution of system of equations (1) and (2) in an iterative manner 
depicted in the sequence of relations (3) to (6). This scheme requires mismatch calculation 
for each half iteration; because [RP] and [RQ] are calculated always using the most recent 
voltage values and it is block Gauss-Seidal approach. The scheme is block successive, 
which imparts increased stability to the solution process. This in turn improves 
convergence and increases the reliability of obtaining solution. 

■ 

PRIOR ART: FAST SUPER DECOUPLED LOADFLOW METHO 
(Refer ences-3, 6, 7) 

Fast Super Decoupled Loadflow (FSDL) Method 

RP P = (APpCosOp + AQj,Sin«&p)/Vp -for PQ-nodes (7) 

* 

RQp = (AQpO>s<Dp-APpSin<I>p)/Vp -for PQ-nodes (8) 
Cosfcp = Absolute (Bpp / /((V + Bj**)) * (" 36 °) (9) 

■ 

Siaflp = -Absolute (Gpp / /(G^ + Bpp 1 )) £ Sin (-36°) (10) 
RP P = AP P / (KpVp) -forPV-nodes (11) 
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PV-nodes or 

(12) 
(13) 



Y9 P p = 2-Ye pq and YVpp = -2V + I-YV M (14) 

bp ,=5 bpCos*p or V = b P (15) 

i 

i 

Kp = Absolute (Bpp/Y9pp) (16) 



Branch admittance magnitude in (12) and (i 3) is of the same algebraic sign as its 
susceptance. Elements of the two gain matrices differ in that diagonal elements of [YV] 
additionally contain the b* values given by relation (15) and in respect of elements 
corresponding to branches connected between two PV-nodes or a PV-node and the slack- 
node. Relations (9) and (10) with inequality sign implies that nodal rotation angles are 
restricted to maximum of -36 degrees. The method consists of relations (3) to (16). In two 
simple variations of the FSDL method, one is to make YV^YBpq and the other is to make 
Y9pq=YVpq. Kp is restricted to the minimum value of 0.75 determined experimentally, and 
it is system independent However it can be tuned for tile best possible convergence for any 
given system. 

This prior art method involves solution of system of equations (1) and (2) in an iterative 
maimer depicted in sequence of relations (3) to (6). Prior art method is embodied in 
algorithm- 1, and in the flow-chart of fig. 1. 

Computation steps of FSDL method (Algorithm-1): 

9 
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Y8pq- 



-Y pq -for branch r/x ratio £ 2.0 

- (Bpq + 0.9(Y pq -Bpq)) -for branch r/x ratio > 2.0 
jBpq -for branches connected between two 

a PV-node and the slack-node 



YVpq = \-Y P q -for branch r/x ratio £ 2.0 

- (Bp q + 0.9(Ypq-Bpq)) -for branch r/x ratio > 2.0 
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■ 

Read system data and assign an initial approximate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

Form nodal admittance matrix, and Initialize iteration count URP = ITRQ= r = 0 
Compute Cosine and Sine of nodal rotation angles using relations (9) and (10), and 
store them. If they, respectively, are less than the Cosine and Sine of -36 degrees, 
equate them, respectively, to those of -36 degrees. 

Form (m+k) x (m+k) size matrices [Y8] and [YV] of (1) and (2) respectively each 
in a compact storage exploiting sparsity. The matrices are formed using relations 
(12), (13), (14), and (15). In [YV] matrix, replace diagonal elements corresponding 
to PV-nodes by very large value (say, 10.0**10). In case [YV] is of dimension (m x 
m), this is not required to be performed. Factorize [Y0] and [YV] using the same 
ordering of nodes regardless of node-types and store diem using the same 
indexing and addressing information* In case [YV] is of dimension (m x m), it is 
factorized using different ordering than that of [Y0]. 

Compute residues [AP] (PQ- and PV-nodes) and [AQ] (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 
Compute the vector of modified residues [RP] using (7) for PQ-nodes, and using 
(11) and (16) for PV-nodes. 

Solve (3) for [A8] and update voltage angles using, [0] = [9] + [AO]. 

Set voltage magnitudes of PV-nodes equal to the specified values, and Increment 

the iteration count ITRP=ITRP+1 and r=(lTRP+ITRQ)/2. 

Compute residues [AP] (PQ- and PV-nodes) and [AQ] (at PQ-nodes only). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 
Compute the vector of modified residues [RQ] using (8) for only PQ-nodes. 
Solve (5) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV]. While 
solving equation (5), skip all die rows and columns corresponding to PV-nodes. 

10 



SUBSTITUTE SHEET (RULE 26) 



WO 2004/023622 PCT/CA2003/001312 

1. Increment the iteration count ITRQ=ITRQ+1 and r=(TrRP+ITRQ)/2, and Proceed 
to step (e) 

m. Calculate line flows and output the desired results 

INVENTED SUPER SUPER DECOUPLED LOADFLOW METHODS 



Super Super Decoupled Loadflow: X'X'-version (SSDLrX'X') 

The general method, in successive (16, IV) iteration scheme represented by sequence of 
relations (3) to (6), can be realized as SSDLOCX 1 , from which manifested are many 
versions. The elements of [RP], [RQ], [Y9] and [YV] are defined by (1 7) to (29). 

RP P = [AP P * + (GppV Bp,,') AQp'] / V p 2 -for PQ-nodes (17) 

RQp = IAQp , -(GppVBpp')APp , ]/Vp -for PQ-nodes (18) 

RPp= [AP P / (Kp*V p 2 )] -forPV-nodes (19) 

■ 

Y6p<i= -1/Xpq* and YV N = -1/Xp,' (20) 



Yepp = E-Ye„ and YVpp = V + I-YVp, (21) 



Where, 

bp* = -2bpCos* p or 
bp 1 = -bpCos* p + [QSHp* - (G w VBppO PSH P '] / V, 2 or 

ty = 2[QSHp> - (Gpp VBppO PSHp'] /V, 2 (22) 

o 

AP P * = APpCos^p + AQpSinOp -for PQ-nodes (23) 

AQ>' = AQpCos$ p - APpSinfy -for PQ-nodes (24) 

ll 
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PSiy = PSHpCosOp + QSHpSin* P -for PQ-nodes (25) 

QSHp* = QSHpCos4> p - PSHpSinfy -for PQ-nodes (26) 

Cos$ p = Absolute [Bpp / v (Gpp 2 + Bpp 2 )] £ Cos (any angle from 0 to -90 degrees) (27) 

♦ 

Sinfcp = -Absolute [Gpp / v (G^, 2 + Bpp 1 )] 2> Sin (any angle from 0 to -90 degrees) (28) 
Kp = Absolute (Bpp/YOpp) (29) 

The factor Kp of (29) is initially restricted to the minimum of 0.75 determined 
experimentally; however its restriction is lowered to the minimum value of 0.6 when its 
average over all PV-nodes is less than 0.6. This factor is system and method independent 
However it can be tuned for the best possible convergence for any given system This 
statement is valid when the factor Kp is applied in the manner of equation (19) in all the 
methods derived in the following from the most general method SSDL-X*X\ 

The definition of Y8p q in (20) is simplified because it does not explicitly state that it always 
takes the value of -Bpq for a branch connected between two PV-nodes or a PV-node and the 
slack-node. This fact should be understood implied in all the definitions of Y8^ in this 
document 



However, a whole new class of methods, corresponding to all those derived in the 
following and prior art, results when the factor Kp is used as a multiplier in the definition of 
RPp at PQ-nodes as in (30) instead of divider in RP P at PV-nodes as given in (19). This 
will cause changes only in relations (17), (19), and (20) as given in (30), (31), and (32). 

RP P - {[AP P ' + (Gpp ' /Bpp') AQp'] / V p 2 } * Kp -for PQ-nodes (30) 

* 

KP P = APp/Vp 2 -for PV-nodes (31) 

12 
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Y9pq - -Bp, and YV pq = -l/X J>q > (32) 

The best performance of methods of this new class has been realized when the factor Kp, 
applied in a manner of relation (30) leading to changes as in (30) to (32), is unrestricted. 
That means it can take any value as given by relation (29). 

■ 

Super Super Decoupled Loadflow: YY-version (SSDL-YY) 

If unrestricted rotation is applied and transformed susceptance is taken as admittance value 
with the same algebraic sign and transformed conductance is assumed zero (reference-6), 
the SSDL-X'X* method reduces to SSDL-YY. Though, this method is not very sensitive to 
the restriction applied to nodal rotation angles, SSDL-YY presented here is the best 
possible experimentally arrived at method However, it gives closely similar performance 
over wide range of restriction applied to the nodal rotation angles (say from -36 to -90 
degrees). 

RP P = AP P V V p 2 and RQp = AQp7V p -forPQ-nodes (33) 

RP P = AP P / (KpVp 2 ) -forPV-nodes (34) 

Ydpq = f-Ypq -for branch r/x ratio £ 3.0 

- (Bp, + (^(Ypq-Bpq)) -for branch r/x ratio > 3.0 
r Bpq -for branches connected between two PV-nodes or 

a PV-node and die slack-node (35) 

YVpq = pYpq -for branch r/x ratio £ 3.0 

\-(Bpq + 0.9(Ypq-Bpq)) -for branch r/x ratio > 3.0 (36) 

Y9pp = Z-YSpq and YVpp - bp' + 2- YV M (37) 

qpT q>p 

bp , = (QSHpVV, l )-b p CosOp or bp* = 2QSHp 7 V,» (38) 

13 
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where, AP P \ AQ P * , QSHp' , and Kp are defined in relations (23) to (29). However, nodal 
rotation angles in relations (27) and (28) be restricted to the maximum of -48 degrees for 
this method, determined experimentally for the best possible convergence from non 
linearity considerations . . 

In case of systems of only PQ-nodes and without any PV-nodes, equations (35) and 
(36) simply be taken as Y6p q = -Ypq and YV pq = -Y M . The factor Kp of (29) is initially 
restricted to the minimum of 0.75 determined experimentally; however its restriction 
is lowered to the minimum value of 0.6 when its average over all PY nodes is less than 
0.6. This factor is system independent However it can be tuned for the best possible 
convergence for any given system 

Branch admittance magnitude in (35) and (36) is of the same algebraic sign as its 
susceptance. Elements of the two gain matrices differ in that diagonal elements of [YV] 
additionally contain the b' values given by relations (37) and (38) and in respect of 
elements corresponding to branches connected between two PV-nodes or a PV-node and 
the slack-node. Relations (27) and (28) with inequality sign implies that nodal rotation 
angles are restricted to maximum of -48 degrees for SSDL-YY. The method consists of 
relation's (3) to (6), (33) to (38), and (23) to (29). In two simple variations of the SSDL- 
YY method, one is to make YVpq-YOpq and the other is to make Y0pq=YVpq. 

SSDL-YY method implements successive (10, IV) iteration scheme and is embodied in 
algorithm-2, and in flow-chart of fig.2 where double lettered steps are characteristic steps 
of the SSDL-YY method and are different than those of the prior art FSDL method 

Computation steps of SSDL-YY method (Algorithm-2): 

a. Read system data and assign an initial approximate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

b. Form nodal admittance matrix, and Initialize iteration count ITRP == ITRQ= r == 0 
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CC 



dd. 



e. 



ff. 



h. 



1. 



J- 



L 



Compute Cosine and Sine of nodal rotation angles using relations (27) and (28), and 
store them. If they, respectively, are less than the Cosine and Sine of -48 degrees, 
equate them, respectively, to those of -48 degrees. 

Form (m+k) x (m+k) size matrices [Y6] and [YV] of (1) and (2) respectively each 
in a compact storage exploiting sparsity. The matrices are formed using relations 
(35), (36), (37), and (38). hi [YV] matrix, replace diagonal elements corresponding 
to PV-nodes by very laige value (say, 10.0**10). In case [YV] is of dimension (m x 
m), this is not required to be performed. Factorize [Y0] and [YV] using the same 
ordering of nodes regardless of node-types and store them using the same 
ing and addressing information. In case [YV] is of dimension (m x m), it is 
factorized using different ordering than that of [Y8]. 

Compute residues AP (PQ- and PV-nodes) and AQ (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 
Compute the vector of modified residues [RP] as in (33) for PQ-nodes, and using 
(34) and (29) for PV-nodes. 

Solve (3) for [AG] and update voltage angles using, [0] « [0] + [A9]. 

Set voltage magnitudes of PV-nodes equal to the specified values, and Increment 

the iteration count ITRP=ITRP+1 andr==(ITRP+ITRQ)/2. 

Compute residues [AP] (PQ- and PV-nodes) and [AQ] (at PQ-nodes only). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 
Compute the vector of modified residues [RQ] as in (33) for only PQ-nodes. 
Solve (5) for [AVJ and update PQ-node magnitudes using [V] = [V] + [AV]. While 
solving equation (5), skip all the rows and columns corresponding to PV-nodes. 
Increment the iteration count ITRQ==ITRQ+1 and i=(ITRP+ITRQ)/2, and Proceed 
to step (e) 

Calculate line flows and output the desired results 



Super Super Decoupled Loadflow: XX-version (SSDL-XX) 

If no (zero) rotation is applied, the SSDL-X'X* method reduces to. SSDL-XX, which is the 
simplest form of SSDL-X'X\ Hie SSDL-XX method comprises relations (3) to (6), (39) to 
(45), and (29). 
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RPp = [AP P + (Gpp / Bpp) AQ P ] / V p 2 -for PQ-nodes (39) 

RQp-tAQp-CCSpp/^APpl/Vp -forPQ-nodes (40) 

RP P = AP P / (KpVp 1 ) -for PV-nodes (4 1) 

Y9pq = [-1 .0/Xpq -for all other branches 

[-Bpq -for branches connected between two PV-nodes or 

a P V-node and the slack-node (42) 



YVpq= -1.0/Xpq -for all branches (43) 



YOpp-E- Y8pq and . W TO = bp' + Z-YV* (44) 

<pp q>p 



bp 1 = -2bp or 

V = -b>+ [QSHp - (Gpp / Bpp) PSHp] /V s 2 or 

V = 2[QSH P - (Gpp / Bpp) PSHp]/V, 2 (45) 



where, Kp is defined in relation (29). This is the simplest method with very good 
performance for distribution networks in absence of PV-nodes (for systems containing only 
PQ-nodes). The large value of the difference [(1/X)-B], particularly for high R/X ratios 
branches connected to PV-nodes, creates modeling error when PV-nodes are present in a 
system. 



Super Super Decoupled Loadflow: BX- version (SSDL-BX) 

If super decoupling is applied only to QV-sub problem, the SSDL-XX method reduces to 
SSDL-BX, which makes it perform better for systems containing PV-nodes. The SSDL- 
BX method comprises relations (3) to (6), (46) to (48), (44) and (45). This method can be 
referred to as Advanced BX-Fast Decoupled Loadflow. 

> 
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RP p - APp/Vp 2 -for aU nodes (46) 

RQp - [AQ P - (Gpp / Bpp) AP P ] / V p -for PQ-nodes (47) 

Y9 pq =-Bp q and YVp^-l/Xp, (48) 



It should be noted that Amerongen's General-purpose Fast Decoupled Loadflow method of 
reference-5 has turned out to be an approximation of this method. The approximation 
involved is only in relation (47). However, numerical performance is found to be only 
slightly better but more reliable than that of the Amerongen's method. 

Super Super Decoupled Loadflow: X'B'-version (SSDLr-X'B') 

RP p = [APp , + (GppVB pp , )AQp > ]/Vp 2 -for PQ-nodes (49) 

RQp = AQpV V p -for PQ-nodes (50) 

RP P = [APp/(Kp*V p 2 )] -forPV-nodes (51) 

Y6pq = -1/Xpq' and YV M = -Bp,* (52) 

YBpp = E-YBpc, and YV„ = V + Z-YVp, (53) 

Where, bp' = -2bpCos<X> p or 

bp , = -bpCos* p + QSHpVV, 3 or 

V = 2QSHpV V. 2 (54) 

Where, AP P \ AQp', PSHp', QSiy, Cosfcp, Sinfcp, K„ are defined in (23)to(29).This method 
consists of relations (3) to (6), (49) to (54), and (23) to (29). Best performance of this 
method could be achieved by restricting <t p in (27) and (28) to less than or equal to -48°. 
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Super Super Decoupled Loadflow: YB'-version (SSDLrYB*) 

The relation (49) in SSOL-XB* implies unrestricted <b p is applied and it can take values up 
to -90 degrees. Therefore, (49) can be modified to (55) with consequent modification of 
(52) into (56). 

RP P = [AP p *Absolute [Bpp / v (G PP 2 + Bpp 1 )] + AQp* [-Absolute [Bpp / v (Gpp 2 + Bpp 1 )]] / V p 2 

-for PQ-nodes (55) 

YGm^-Yp, and YV pq = -B M ' (56) 

This method consists of relations (3) to (6), (55), (50), (51), (56), (53) and (54), and (23) to 
(29). Best performance of this method could be achieved by restricting 3> p in (27) and (28) 
to less than or equal to -48 degrees. Where, Aiy , AQ P ' , PSHp' , QSHp' , Cos$ p , Sin* p , 
Kp are defined in (23) to (29). 

Super Super Decoupled Loadflow: B'X'-version (SSDL-B'X') 

RP P - AP P 7 V p 2 -for PQ-nodes (57) 

RQ„ = [AQp' - (GppV B„,') AP P '] / V p -for PQ-nodes (58) 

RP P = [APp/^Vp 2 )] -forPV-nodes (59) 

Y %i = -Bp,' and W^-lfX^ (60) 

Yepp = Z-Y9p, and YVpp = V + Z-YV M (61) 

<f*p <pp 

Where, bp' = -2bpCos3> p or 

V = -bpCos4>p+[QSHp , -(GppVB P p')PSHp > yV. 2 or 
V=2[QSHp'-(C3j,VBpp')PSHp']A^. 2 (62) 
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Where, AP P \ A<V, PSHp\ QSiy, Cos* p , SinO p , Kp are defined in (23)to (29).This 
method consists of relations (3) to (6), (57) to (62), and (23) to (29). Best performance of 
this method could be achieved by restricting <D P in (27) and (28) to less than equal to -48°. 

Super Super Decoupled Loadflow: B'Y-version (SSDL- B'Y) 

The realtion (58) in SSDL-B'X' implies unrestricted O p is applied and it can take values up 
to -90 degrees. Therefore, (58) can be modified to (63) with consequent modification of 
(60) into (64). 

RQp - [AQp*Absolute [Bpp / v (Gpp 2 + B^)] - AP p * [-Absolute [Bpp I v ((V + Bpp 2 )]] / V p 2 

-for PQ-nodes (63) 
YOpq - -Bp,' and YV M = -Y w (64) 

This method consists of relations (3) to (6), (57), (63), (59), (64), (61) and (62), and (23) to 
(29). Best performance of this method could be achieved by restricting <D P in (27) and (28) 
to less than or equal to -48 degrees. Where, AP P ' , AQ P ' , PSiy , QSiy , CosO p , SinO p , 
Kp are defined in (23) to (29). 

Simultaneous (IV, 1 8) Iteration Scheme 

An ideal to be approached for the decoupled Loadflow methods is the constant matrix 
Loadflow of reference-6 referred in this document as BGGB -method In an attempt to 
imitate it, a decoupled class of methods with simultaneous (IV, 16) iteration scheme 
depicted by sequence of relations (65) to (69) is invented This scheme involves only one 
mismatch calculation in an iteration. The correction vector is calculated in two separate 
parts without intermediate updating. Each iteration involves one calculation of [RQ], [AVJ, 
and [RP], [ A9] to update [V] and [6]. 

[AV] = [YVy 1 [RQ] (65) 
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[RP] = [AP/V] - [G] [AV] (66) 



[A9] - [Y9] 1 [RP] (67) 



[0] = [6] + [A8] (68) 



[V] = [V] + [AV] (69) 



In this invented class, each method differs only in the definition of elements of [RQ] and 
[YV]. The accuracy of methods depends only on the accuracy of calculation of [AV]. The 
greater the angular spread of branches terminating at PQ-nodes, the greater the inaccuracy 
in the calculation of [ AV] . 

Super Super Decoupled Loadflow: BGX'-version (SSDLrBGX*) 

Numerical performance could further be improved by organizing the solution in a 
simultaneous (IV, 10) iteration scheme represented by sequence of relations (65) to (69). 
The elements of [RP], [RQ], [Y0] and [YV] are defined by (70) to (74). 

RQp = [AQp , -(Gpp , /Bpp')AP p , ]/Vp -for PQ-nodes (70) 

RP P = (APp/Vp^ZGpqAVq -for all nodes (71) 

YSp, = -Bp, and YVp, = -1/Xp,* (72) 



Y9pp = i:-Yep q and YVp, = V + E-YVp, (73) 



bp' = -2bpCos* p or 
bp^-tpCosOp + IQSHp'^G^VBpp'JPSI^'yV, 1 or 

V = 2[QSHp' - (GppV Bpp') PSHp'V V. 2 (74) 
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Where, AP P ' , AQp* , PSHp' , QSH P ' , Cos<D p , Sin<D p are defined in (23) to (28).The SSDL- 
BGX' method comprises relations (65) to (74), and (23) to (28). Best possible convergence 
could be achieved by restricting rotations (<& p ) in the range (-10° to -20°) in relations (27) 
and (28). The method is embodied in algorithm-3 and in the flow-chart of Fig.3. 

Super Super Decoupled Loadflow: BGY- version (SSDL-BGY) 

If unrestricted rotation is applied and transformed susceptance is taken as admittance 
values and transformed conductance is assumed zero (reference-6), the SSDI^BGX* 
method reduces to SSDL-BGY as defined by relations (75) to (79). 

RQ P = AQpWp = (AQ p Cos<I> p - AP p Sin0> p )/V p -for PQ-nodes (75) 

m 

RP P = (APj/Vp) - ZGpqAV, -for all nodes (76) 

Yep, = -Bp, and YV M = -Y pq (77) 



Y6pp = Z-YBp, and YVpp = V + S-YV„ 

<F*P q>p 

bp = -2bpCos* p or 
bp' = -bpCos*p + (QSHpCos0 p - PSHpSin*p)/V. 2 or 
V = 2(QSHpCos* p - PSHpSinO p yV, 2 

♦ 

The SSDL-BGY method comprises relations (65) to (69), and (75) to 
case of the SSDL-BGX' method 

Super Super Decoupled Loadflow: BGX-version (SSDL-BGX) 

♦ 

If no (zero) rotation is applied, the SSDL-BGX* method reduces to SSDL-BGX as defined 
by relations (80) to (84). 
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RQp = [AQp - (Gpp / Bpp)AP P 3 / V p -for PQ-nodes 



zn 



RP P = (APpA^ p )-EGp,AV q -for all nodes 



Y0p q - -Bp, 



and 



YVpq -1/Xpq 
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(80) 



(81) 



(82) 



Y8pp = S-YSp, 



and 



YVpp = bp' + Z-YVp, 

q>p 



bp' = -2bpCosO p 

V = -bpCos* p + [QSHp - (Gpp/Bpp)PSH,yV, 2 

V = 2[QSHp - (Gpp/Bpp)PSHpyV. 2 



or 



or 



(83) 



(84) 



The SSDL-BGX method comprises relations (65) to (69), and (80) to (84). It is again the 
special case of the SSDL-BGX' method. 

Computation steps of SSDL-BGX', SSDL-BGY and SSDL-BGX methods 

* 

(Algorithm-3): 



b. 

ccc 



Read system data and assign an initial approximate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 
Form nodal admittance matrix, and Initialize iteration count UR = 0. 
Compute Sine and Cosine of nodal rotation angles using relations (28) and (27), and 
store them. If they, respectively, are less than the Sine and Cosine of any angle set 
(say in the range -10 to -20 degrees), equate than, respectively, to those of the same 
angle in the range -10 to -20 degrees. In case of zero rotation, Sine and Cosine 
value vectors are not required to be stored 
dd<L Form (m+k) x (nrfk) size matrices [Y6] and [YV] of (l) and (2) respectively each 
in a compact storage exploiting sparsity 

1) In case of SSDL-BGX'-method, the matrices are fonned using relations (72), 
(73), and (74) 
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2) In case of SSDL-BGY-method, the matrices are formed using relations (77), 
(78), and (79) 

3) In case of SSDL-B GX-method, the matrices are formed using relations (82), 
(83), and (84) 

In [YV] matrix, replace diagonal elements corresponding to PV-nodes by very large 
value (say, 10.0**10). In case [YV] is of dimension (m x m), this is not required to 
be performed. Factorize [Y0] and [YV] using the same ordering of nodes regardless 
of node-types and store them using the same indexing and addressing 
information. In case [YV] is of dimension (m x m), it is factorized using different 
ordering than that of [Y9]. 

e. Compute residues AP (PQ- and PV-nodes) and AQ (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

£f£ Compute the vector of modified residues [RQ] using (70) in case of SSDL-BGX? , 
using (75) in case of SSDL-BGY, and using (80) in case of SSDL-BGX method for 
only PQ-nodes. Solve (65) for [AV]. While solving equation (65), skip all the rows 
and columns corresponding to PV-nodes. Compute the vector of modified residues 
[RP] using (71) or (76) or (81). Solve (67) for [AG], 

ggg. Update voltage angles using, [0] = [9] + [A8]. and update PQ-node voltage 
magnitudes using [V] = [V] + [AV]. 

hhh. Set voltage magnitudes of PV-nodes equal to the specified values, and Increment 
the iteration count ITR=ITR+1 , and proceed to step (e). 

del Calculate line flows and output the desired results 

Triple lettered steps are characteristic steps of algorithm-3. The SSDL-BGX*, SSDL-BGY 
and SSDL-BGX methods differ only in steps-ccc and -ddd defining gain matrices, and 
step-flf for calculating [RP] and [RQ]. Fig. 3 is the flow-chart embodiment of algorithm-3. 

Super Super Decoupled Loadflow: X'G^X'-version (SSDL- X'GpyX 9 ) 

Numerical performance could also be improved by organizing the solution in a 
simultaneous (IV, 19) iteration scheme represented by sequence of relations (65) to (69). 
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The elements of [RP], [RQJ, [Y6] and [YV] for this method are defined by (85) to (91). 



RP P = {[AP P ' + (Gpp 7 BppOAQ,'] / V p 2 } - (gp'AV p ) -for PQ-nodes (85) 



RQ p = [AQ p , -(G P p/Bpp , )APp , ]/Vp -for PQ-nodes (86) 

RP P = [(APpA^p 2 ) - £ Gp, AVJ / Kp -forPV-nodes (87) 

Y6 pq = -1/Xp q ' and YVp^-l/Xp,' (88) 



Ye pp = Z-Ye pq and YV PP = bp' + S-YVp, (89) 

<pp q>p 



bp' = -2bpCos* p or 
V = -bpCosOp + [QSHp* - (GppV BppOPSHp']^. 2 or 

bp' = 21QSH,' - (Gpp 7Bpp')PSHp']/V, 2 (90) 

♦ 

gp' = 2b^Sin* p or 
gp" = bpSin<&p + [PSHp' + (GppV BppOQSHp'yVs 2 or 

gp' - 2[PSJV + (Gpp 7 B^OQSIVJ/V, 2 (91) 

Where, AP P * , AQp* • PSiy , QSiy , Cos* p , Sin* p , Kp are defined in (23) to (29) Again, 
if unrestricted rotation is applied and transformed susceptahce is taken as admittance values 
and transformed conductance is assumed zero (reference-6), the SSDL-X'G^X' method 
reduces to SSDL-YG^Y. If no (zero) rotation is applied, the SSDLrX'GpvX' method 
reduces to SSDL-XGpvX. The SSDL-X'G^X' method comprises relations (65) to (69), 
(85) to (91), and (23) to (29). It is embodied in algorithm-4 and in the flow-chart of Fig.4. 

Computation steps of SSDL-X'G^X', SSDL-YG^Y and SSDL-XG^X methods 
(Algorithm-4): 
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a. Read system data and assign an initial approximate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

b. Form nodal admittance matrix, and Initialize iteration count ITR= 0. 

ccc Compute Sine and Cosine of nodal rotation angles using relations (28) and (27), 
store them. If they, respectively, are less than the Sine and Cosine of any angle set 
(say 0 to -90 degrees), equate them, respectively, to those of the same angle in the 
range 0 to -90 degrees. In case of zero rotation, Sine and Cosine vectors are not 
required to be stored. 

dddd. Form (m+k) x (m+k) size matrices [Y0] and [YV] of (I) and (2) respectively each 
in a compact storage exploiting sparsity using relations (88), (89), and (90). 
In [YV] matrix, replace diagonal elements corresponding to PV-nodes by very large 
value (say, 10.0**10). hi case [YV] is of dimension (m x m), this is not required to 
be performed. Factorize [Y8] and [YV] using the same ordering of nodes regardless 
of node-types and store them using the same indexing and addressing 
information. In case [YV] is of dimension (m x m), it is factorized using different 
ordering than that of [Y8]. 

e. Compute residues AP (PQ- and PV-nodes) and AQ (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

ffiC Compute [RQ] using (86) for only PQ-nodes. Solve (65) for [AV]. While solving 
equation (65), skip all the rows and columns corresponding to PV-nodes. Compute 
the vector of modified residues [RP] using relations (85), (87), and (29). Solve (67) 
for [AG]. 

ggg. Update voltage angles using, [6] = [9] + [A6]. and update PQ-node voltage 

magnitudes using [Vj = [V] + [AV], 
hhh. Set voltage magnitudes of PV-nodes equal to the specified values, and Increment 

the iteration count ITR=ITR+1 , and proceed to step (e) 
m. Calculate line flows and output the desired results 
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Four lettered steps are characteristic steps of algorithm-4. This method is useful particularly 
for distribution systems without PV-nodes. Fig.4 is the flow-chart embodiment of 
algorithm-4. 



Common Statements Concerning all methods: 

In ail the prior art and invented models [Y9] and [YV] are real, sparse, symmetrical and 
built only from network elements. Since they are constant, they need to be factorized once 
only at the start of the solution Equations (1) and (2) are to be solved repeatedly by 
forward and backward substitutions. 

[Y8] and [YV] are of the same dimensions (m+k) x (m+k) when only a row/column of the 
slack-node is excluded and both are triangularized using the same ordering regardless of 
the node-types. For a row/column corresponding to a PV-node excluded in [YV], use a 
large diagonal to mask out the effects of the off-diagonal terms. When the node is switched 
to the PQ-type the row/column is reactivated by removing the large diagonal. This 
technique is especially useful in the treatment of PV-nodes in the matrix [YV] . 

It is invented to make this technique efficient while solving (5) or (65) for [AV] by 
skipping all PV-nodes and fector elements with indices corresponding to PV-nodes. In 
other words efficiency can be realized by skipping operations on rows/columns 
corresponding to PV-nodes in the forward-backward solution of (5) or (65). Ihis has been 
implemented and time saving of about 4% of the total solution time (including 
input/output) could be realized in 14-14 iterations required to solve 118-node system with 
the uniform R-scale factor of 4 applied It should be noted that the same indexing and 
addressing information can be used for (he storage of both matrices as they are of the 
same dimension and spar si ty structure. 

ALGORITHMS using GLOBAL CORRECTIONS 

The algorithms- 1, -2, -3, and -4 in the above involve incremental (or local) ccarections. All 
the above algorithms can be organized to produce corrections to the initial estimate 
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solution. It involves storage of the vectors of modified residues and replacing the relations 
(17), (18), (19) by (92), (93), (94) respectively, and (4) or (68) and (6) or (69) respectively 
by (95) and (96). Superscript '0 1 in relations (95) and (96) indicates the initial solution 



estimate. 

RV - KAPpT + (Gpp 7 BppO (AQp 1 ) ']/(\V) 2 + RP P (M > (92) 

RQp r - [(AQp 1 )* - (GppVIV) (APpO'VOV) 2 + RQp (M) (93) 

■ 

RP P r = AV/tKpOVfi + RPpfr- 9 (94) 

ep r =e p 0 +Aep r (95) 

V p r = V P ° + AVp r (96) 



RECTANGULAR COORDINATE FORMULATIONS OF INVENTED 
LOADFLOW METHODS 

1 

This involves following changes in the equations describing the loadflow models 
formulated in polar coordinates. 

(i) Replace 6 and A9 respectively by f and Af in equations (1),(3), (4), (67), (68) and (95) 

(ii) Replace V and AV respectively by e and Ae in equations (2), (5), (6), (65), (66), (69) 
and (96) 

(iii) Replace V p by e, and V, by e, in equations (17) to (19), (22), (30), (31), (33), (34), 
(38) to (41), (45) to (47), (49) to (51), (54), (55), (57) to (59), (62), (63), (70), (71), 
(74) to (76), (79) to (81), (84) to (87), (90), (91). The subscript V indicates the slack- 
node variable. 
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(iv) After calculation of corrections to the imaginary part of complex voltage (Af) of 
PV-nodes and updating the imaginary component (f) of PV-nodes, calculate real 
component by: 



ep =>|V-f p 2 



(97) 



APPENDIX 



Transformation of Branch Admittance 

The branch admittance transformation for symmetrical gain matrices of the methods 
described in the above is given by the following steps: 

1 . Compute: <X> P = arctan (Gpp/Bpp) and 

<& q = arctan (G qq /B qq ) (98) 

2. Compute the average of rotations at the terminal nodes (p and q) of a branch: 

* 

Oar^^p + O)^ (99) 

3. Compare with the Limiting Rotation Angle (LRA) and let * a v to be the smaller 
of the two: 



O av = minimum (* 8 *, LRA) (i oo) 



Compute transformed pq-th element of the admittance matrix* 



<V + J*W " (COS + j Sin O av ) (Gpq + jBpq) (101) 
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5. Note that the transformed branch reactance is: 

» a 

Xpq^Bpq^Gpq^ + Bpq* 2 ) and similarly, (102) 

:V=B PP y(Gp 2 +B PJ) ' 2 ) (103) 

In the description above Xpq* is the transformed branch reactance defined by equation (1 03) 
and Bpq* is the corresponding transformed element of the susceptance matrix Gpp 1 and Bpp* 
are diagonal elements obtained from (102). 

SOME POSSIBLE SIMPLE VARIATIONS OF SSDL-METHODS 

1. Simple obvious modifications are the use of V p and V p 2 interchangeably in all 
expressions of RP P , and the use of 1 .0 for V 9 2 in all expressions of bp 9 involving the 
termV, 2 

2. can also take values without transformation of bp and QSH P 

3. Explicit algorithmic steps are not given for many variants of SSDL-X'X* except 
SSDL-YY, They are obvious from their descriptions and are similar to those of 
SSDL-YY method 

EXPLANATORY STATEMENTS 

The system stores a representation of the reactive capability characteristic of each machine 
and these characteristics act as constraints on the reactive power, which can be calculated 
for each machine. 

While the description above refers to particular embodiments of the present invention, it 
will be understood that many modifications may be made without departing from the spirit 
thereof. The accompanying claims are intended to cover such modifications as would fall 
within the true scope and spirit of the present invention. 
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The presently disclosed embodiments are therefore to be considered in all respect as 
illustrative and not restrictive, the scope of the invention being indicated by the appended 
claims in addition to the foregoing description, and all changes which come within the 
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therein. 
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The present invention is applicable to systems to process Loadflow computation by means 
of modified real and reactive power residues, and gain matrices derived from the Jacobian 
matrix. Hie embodiments of the invention in which an exclusive property or privilege is 
claimed are defined as follows: 

1 . A method of controlling security (over load, under/over voltage) in a power system, 
comprising the steps of: 

obtaining on-line data of nodal injections, voltages and phases at main nodes, and 
open/close status of circuit breakers in the power system, 

establishing initial specifications of controlled parameters (real and reactive power 
at PQ-nodes, real power and voltage magnitude at PV-nqdes, and 
transformer turns ratios etc.), 

performing Loadflow computation at said nodes of the power system by a Super 
Super Decoupled computation in any of the Super Super Decoupled 
Loadflow methods or any of their hybrid combination or simple variants 
employing corresponding gain matrices derived from a super decoupled 
Jacobian matrix for real power with respect to angle and a super decoupled 
Jacobian matrix for reactive power with respect to voltage, and involving 
triangular factorization of said gain matrices and computing of discrepancy 
of real powa* and reactive power from specified values through such nodes, 

said computing including introducing variables representing quotients of the 
transformed discrepancies from specified values of real and reactive power 
flowing in through each node divided by voltage, or square of the voltage in 
case of transformed real power mismatches in methods employing (19, IV) 
iteration scheme, on each node, and using such variables to calculate values 
of angle and voltage for said transformed discrepancies from specified 
values of real and reactive power flowing in thro ugh each node, by using 
triangular factorization of said gain matrices for real and reactive power, 
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» 

initiating said Loadflow computation with guess solution of the same voltage 
magnitude and angle as those of the slack (reference) node referred to as 
slack start, 

restricting nodal transformation angle to maximum -48 degrees, applied to complex 
power injection in computing said transformed discrepancies from specified 
values of real and reactive power flowing in through each node, 

evaluating the computed Loadflow for security (over load, under/over voltage), 

correcting one or more controlled parameters with said correction (amount of over 
load and/or undo/over voltage) values and repeating the computing and 
evaluating steps until evaluating step finds a good power system (no over 
load, no under/over voltage), and 

effecting a change in the voltages and phases at said nodes of the power system by 
actually implementing the finally obtained values of controlled parameters 
after evaluating step finds a good power system. 

A method as defined in claim 1 wherein said Super Super Decoupled methods, 
employing successive (18, IV) iteration scheme, of Loadflow computation 
are characterized in modifying the transformed real power residue at a PQ- 
node and real power residue at a PV-node by dividing them by squared 
voltage magnitude multiplied by a factor determined as a ratio of a diagonal 

» - * 

element of susceptance matrix to a diagonal element of corresponding said 
gain matrix derived from transformed Jacobian matrix for real power with 
respect to angle. 

A method as defined in claiml wherein said Super Super Decoupled methods, 
employing simultaneous (IV, 10) iteration scheme, of Loadflow 
computation are characterized in that they involve only one time calculation 
of real and reactive power residues in an iteration unlike two calculations in 
successive (16, IV) iteration scheme, and thereby achieving consequent 
speed-up. 

A simple system/method of controlling generator and transformer voltages of more 
elaborate method of security control defined in claim 1 can be realised in a 
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system for controlling generator and transformer voltages in an electrical 
power utility containing plurality of electromechanical rotating machines, 
transformers and electrical loads connected in a network, each machine 
having a reactive power characteristic and excitation element which is 
controllable for adjusting the reactive power generated or absorbed by the 
machine, and some of the transformers having controllable taps for adjusting 
terminal voltage, said system comprising: 

means de fining any of Super Super Decoupled models of the network for providing 
an indication of the quantity of reactive power to be supplied by generators 
including at a reference node in dependence on representations of selected 
network electrical parameters, 

machine control means connected to the said excitation element of at least one of the 
rotating machines for controlling the operation of the excitation element of 
at least one machine to produce or absorb the amount of reactive power 
indicated by said model means with respect to the set of representations. 

5. A system as defined in claim 4 wherein the network includes a plurality of nodes 

each connected to at least one of: a reference generator, a rotating machine; 
and an electrical load, and the model has one of the 3 -forms and their 
variants of Super Super Decoupled matrices which receives representations 
of selected values of the real and reactive power flow from each machine 
and to each load, and the model is operative for producing calculated values 
for the reactive power quantity to be produced or absorbed by each machine. 

6. A system as defined in claim 4 wherein the utility further has at least one transformer 

having an adjustable transformation ratio, and said means defining a model 
is further operative for producing a calculated value for the transformer 
transformation ratio. 

7. A system as defined in claim 4 wherein said machine control means are connected to 

said excitation element of each machine for controlling the operation of the 
excitation element of each machine. 
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A system as defined in claim 4 wherein said transformation ratio control means are 
connected to said transformer tap changing element of each transformer for 
controlling the operation of the transformer tap changing element of each 
transformer. 

A method far controlling generator and transformer voltages in an electrical power 
utility containing plurality of electromechanical rotating machines, 
transformers and electrical loads connected in a network, each machine 
having a reactive power characteristic and excitation element which is 
controllable for adjusting the reactive power generated or absorbed by the 
machine, and some of the transformers having controllable taps for adjusting 
terminal voltage, said method comprising: 

creating any of said Super Super Decoupled models of the network for providing an 
indication of the quantity of reactive power to be supplied by the, generators 
in dependence on representations of selected network electrical parameters, 

controlling the operation of the excitation element of at least one machine to 
produce or absorb the amount of reactive power indicated by any of the said 
Super Super decoupled models with respect to the set of specified 
parameters. 

A method as defined in claim 9 wherein said step of controlling is carried out to 

control the excitation element of each machine. 
A method as defined in claim 9 wherein said step of controlling is carried out to 

control the tap-changing element of each transformer. 
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